Investigation of mechanical properties of twin gold crystal nanowires under uniaxial load by molecular dynamics method
Zhang Guo-Wei, Yang Zai-Lin†, , Luo Gang
College of Aerospace and Civil Engineering, Harbin Engineering University, Harbin 150001, China

 

† Corresponding author. E-mail: yangzailin00@163.com

Project supported by the National Science and Technology Pillar Program, China (Grant No. 2015BAK17B06), the Earthquake Industry Special Science Research Foundation Project, China (Grant No. 201508026-02), the Natural Science Foundation of Heilongjiang Province, China (Grant No. A201310), and the Scientific Research Starting Foundation for Post Doctorate of Heilongjiang Province, China (Grant No. LBHQ13040).

Abstract
Abstract

Twin gold crystal nanowires, whose loading direction is parallel to the twin boundary orientation, are simulated. We calculate the nanowires under tensile or compressive loads, different length nanowires, and different twin boundary nanowires respectively. The Young modulus of nanowires under compressive load is about twice that under tensile load. The compressive properties of twin gold nanowires are superior to their tensile properties. For different length nanowires, there is a critical value of length with respect to the mechanical properties. When the length of nanowire is greater than the critical value, its mechanical properties are sensitive to length. The twin boundary spacing hardly affects the mechanical properties.

1. Introduction

In recent years, metal nanomaterials, whose properties are different from macroscopic metal materials, have received a great deal of attention. The experiment of metal nanowires is uncontrollable, but by the molecular dynamics (MD) method their deformation can be observed under a certain condition.

Gleiter[1,2] investigated the deformation of nano-structure. Janusz et al.[3] used the molecular dynamics simulation method to calculate thermodynamic properties of the transition metal with face-centered cubic structure. Park and Zimmerman[4] discussed the mechanical properties of the gold nanowires in the yield condition, and summarized the plastic characteristics of gold nanowires by molecular dynamics simulation.

The current researches focus on the twin structure, whose load direction is vertical to the twin boundary orientation. The point defects in Au twin boundaries,[5] the fracture behavior and the ductility of the twinned Cu nanowires,[6] and the ductility of the magnesium nanowire[7] were studied under loads in the above direction.

In the present paper, we discuss the properties of gold twin crystal nanowires, whose directions are all parallel to the twin boundary orientation, under different conditions.

The rest of the paper is organized as follows. In Section 2, we discuss the simulation model. The results will be discussed in Section 3. Finally, some conclusions are drawn from the present study in Section 4.

2. Simulation model

Some twin gold crystal models are given by the geometric construction method (Fig. 1). The nanowire models are rectangular columnar, whose xy cross-section is the twin crystal section (Fig. 1(a)). The transverse section sizes of models are all 8a × 8a, of which a is the gold lattice constant and a = 0.4078 nm = 4.078× 10−10 m. The other two cross-sections are shown in Figs. 1(b) and 1(c). The lengths of nanowires are all 19a and the atoms at both ends of the z axis are tool atoms, which are under loads and not involved in statistical computation.

Fig. 1. Three view drawings of the simulation model.

The molecular dynamics simulation software NanoMD,[8] which is developed by Zhao’s group in Nanjing University, are modified[9,10] for simulation. In this paper we adopt Johnson’s EAM potential function[11] to describe the interaction between each two gold atoms. Cell list and Verlet list are combined with integration step δt = 2.94 fs = 2.94× 10−15 s. At the beginning, atom initial velocity, which fits the Maxwell–Boltzmann distribution,[12] is generated randomly. The temperature of the model system keeps up 300 K by the Nosé-Hoover[1316] method. Through 2×104 relaxation steps, the system achieves a thermal equilibrium state. The loads is exerted on the [110] orientation in the z-axial direction, whose stress is obtained by[17]

where is the zz component of the atomic stress tensor of atom α; Ωα is its volume; mα is the mass; is the velocity component in the z direction of atom α; rα β is the distance between atom α and atom β; ϕ, F, ρ, and f are parameters from EAM potential.

3. Result and discussion

The nanowires under tensile or compressive loads and different lengths are simulated respectively. Their mechanical properties are different under different conditions.

3.1. Properties of twin crystal nanowires under different load styles

The models deform at 0.129 m/s tensile or compressive velocity rate in the z-axial direction on the basis of the above mentioned models.

3.1.1. Energy discussion

The average energies of models under different load styles are different (Fig. 2). The average kinetic energy (KE) of the system is as follows:

where kB, N, T are the Boltzmann constant, the number of atoms, and the temperature of the system respectively. When the system is definite, the average kinetic energy is only depending on the temperature of the system. The kinetic energies of models are almost the same, whose error is less than 1× 10−5 under tensile or compressive load (Fig. 2(a)).

Fig. 2. Kinetic energy (a) and potential energy (b) of nanowires under tensile and compressive loads.

The average potential energies in different steps under tensile and compressive loads are shown in Fig. 2. With ignoring the accidental error of the computer simulation, both relaxation processes are the same. After relaxation, the potential energy increases with load. Under tensile load the potential energy of the system is almost constant until the 105-th step, however under compressive load the potential energy increases in an oscillatory manner.

3.1.2. Mechanical properties discussion

Figure 3 shows that with the strain increasing each of the curves is almost a line like that in the elastic stage of macroscopic material at the beginning of deformation. Then the system arrives at the yield stage and the stress at the peak of the curve is defined as the yield stress σs. After yielding, system stress decreases rapidly and its speed is faster under tensile load than under compressive load.

Fig. 3. Relationship between stress and strain of nanowires under compressive load and tensile load.

Table 1 gives some mechanical values of twin gold crystal nanowires under two load styles, where σs, ɛs, E are yield stress, yield strain and Young modulus respectively. The nanowires yield more easily under compressive load than under tensile load, however the Young modulus in the former case is twice as great as that in the latter case.

Table 1.

Mechanical values of nanowires under different load styles.

.
3.1.3. Micro-deformation discussion

In the paper, we discuss the micro-deformations by using colored atom balls, where blue atoms are the twin boundary or dislocation atoms and yellow atoms are stable.

Relaxations under tensile and compressive loads are similar due to the fact that models are not applied loads. In the process of compression, compressive instability leads to the fact that dislocations appear, which come into slip planes. The models will be broken in the form of pileup along slip planes. Nevertheless, there are different statuses in the models under tensile loads. Slip planes break the twin crystal structure, around which the non-crystallizing of the atoms is acute and the necking phenomenon turns up until it is broken. The directions of dislocations are both 45 angle and twin boundaries prevent the dislocation from developing.

Fig. 4. Micro-deformations of nanowires under compressive load and tensile load.
3.2. Properties of different length nanowires

On the basis of the above mentioned models, different length twin gold nanowires are simulated. Different mechanical properties are obtained.

3.2.1. Energy discussion

The energy increases linearly during the relaxation. The system’s average kinetic energy is dependent only on temperature (Eq. (2)), so the average kinetic energies of different length models are coinciding with each other in the process of the whole simulation. The potential energy curves are different while their tendencies are similar. Also, the longer the nanowire, the lower the potential energy curve is (Fig. 5(a)). Figure 5(b) shows the potential energy at a certain step, where Ef is the potential energy at the final step of the simulation. All certain potential energies are small with the length increasing, and the final step potentials are biggest, which verifies the fact that the potential energy of the steady state is smallest.

Fig. 5. (a) Average potential energies versus simulation step and (b) the particular step potential energies versus simulation step for different length nanowires.
3.2.2. Mechanical properties discussion

Although the lengths of gold nanowires are different, the variation trends of curves are similar (Fig. 6). The system begins with elastic deformation until the yield stage, then falls in an oscillatory manner until the break. The longer the nanowire, the more rapidly the break happens. Table 2 gives the mechanical values of different length twin gold crystal nanowires in the yield stage. When the length L > 19a, Young modulus E diminishes clearly and the curves (Fig. 6) have no secondary peak values.

Table 2.

Mechanical values of different length nanowires in the yield step.

.
Fig. 6. Stress–strain curves of different length nanowires.
3.2.3. Micro-deformation discussion

Figure 7 shows representative microscopic structures of different length twin gold nanowires in broken stages. When the nanowires yield, low coordination number atoms of the system increase. Then slip planes come out immediately, at this point, stress decreases with the increase of strain. Twin gold nanowires with a shorter length form a stress well, because most of the atoms in the nanowires are in a non-crystallizing state and there are many low coordination number atoms. Through a combination of low coordination atoms and slip planes, a stress well appears in each of the stress–strain curves of the nanowires. As the simulation continues, the non-crystallizing atomic motion of the atoms around slip planes increases, which brings about the damage of the initial twin crystal structure. Then the necking phenomenon appears, resulting in the fracture of nanowires. Under the same tensile speed load, the longer the nanowires, the smaller the influence of non-crystallizing atomic motion of the atoms around slip planes on the whole system after yielding.

Also, in the progress of tension it is verified that the dislocations destroy the prevention of a twin boundary and keep developing along the length direction.[18]

Fig. 7. Micro-deformations of different length nanowires.
3.3. Properties of different twin boundary nanowires

The twin boundary spacing has an effect on the mechanical properties when the tensile direction is vertical to the twin boundary orientation.[19] There appear different properties when the loading direction is parallel. The other two twin boundary models are simulated (Fig. 8). Because only twin boundary spacings are different, the average energies are all the same.

Fig. 8. Three twin boundary nanowires (a) and their stress–strain curves (b).

Figure 8(b) shows that the yield stresses of nanowires under tensile loading, whose directions are parallel to the twin boundary, are the same, which is different from the scenario under vertical tensile loading. Only the stress of the 5 twin nanowire decreases faster after yielding. The twin boundary spacing affects only the rate of dislocation developing.

4. Conclusions

In this work, we investigate the mechanical properties of twin gold crystal nanowires by the molecular dynamics method through controlling variables. The nanowires under tensile or compressive loads, different length nanowires, and different twin boundary nanowires are simulated respectively. The micro-deformation is the progress of dislocation developing by breaking twin boundaries.

When the models are under different style loads, the mechanical properties are different. The Yong modulus under a compressive load is about twice that under a tensile load. The compressive properties of twin gold nanowires are superior to its tensile properties.

The length of nanowires also affects the mechanical properties. When the length of nanowires is short: L = 10a, 13a, 16a, the nanowires length changes hardly affect the mechanical properties of twin gold nanowires; when the length turns long with the increase of the nanowire length, system yield stress and Young modulus decrease and the yield strain increases. The longer the nanowires, the less the influence of non-crystallization is.

When the loading direction is parallel to the twin boundary, the twin boundary spacing hardly affects the mechanical properties.

Reference
1Gleiter H 2000 Acta Mater. 28 1
2Gleiter H 1995 Nanostruct. Mater. 6 3
3Diao J KGall KDunn M L2003Nat. Mater.552465
4Park H SZimmerman J A 2005 Phys. Rev. 72 054106
5Hussain FHayat S SShah Z AHassan NAhmad S A 2013 Chin. Phys. 22 096102
6Sun JFang LMa AJiang JHan YChen HHan J 2015 Mater. Sci. Eng. A, Struct. Mater.: Prop. Microstruct. Process. 634 86
7Ni CDing HJin X 2016 Comp. Mater. Sci. 111 163
8Zhao JYin XLiang SLiu YWang DDend SHou J 2008 Chem. Res. Chin. Univ. 24 367
9Yang ZZhang GLuo G 2015 Chin. Phys. 24 066203
10Yang ZZhang GZhao J 2016 Phys. Lett. 380 917
11Johnson R A 1989 Phys. Rev. 39 12554
12Long ZXu WTang L2007Molecular Simulation Theory and PracticeVol. 1BeijingChemical Industry Press798179–81(in Chinese)
13Hoover W G 1985 Phys. Rev. 31 1695
14Nośe S 1984 Chem. Phys. 81 511
15Nośe S 1984 Mol. Phys. 52 255
16Hoover W G 1985 Phys. Rev. 31 1695
17Wu H A 2006 Eur. J. Mech. A, Solids 25 370
18Cao AWei Y 2006 Phys. Rev. 74 214108
19Cao A JWei Y GMao S X 2007 Appl. Phys. Lett. 90 151909